Neutron Scattering Investigation of the Spin Ice State in Dy2Ti207. 
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Dy2Ti207 has been advanced as an ideal spin ice material. We present a neutron scattering 
investigation of a single crystal sample of ^^^Dy2Ti207. The scattering intensity has been mapped 
in zero applied field in the h, h, I and h, k, planes of reciprocal space at temperatures between 
0.05 K and 20 K. The measured diffuse scattering has been compared with that predicted by the 
dipolar spin ice model. The comparison is good, except at the Brillouin zone boundaries where 
extra scattering appears in the experimental data, ft is concluded that the dipolar spin ice model 
provides a successful basis for understanding Dy2Ti207, but that there are issues which remain to 
be clarified. 



PACS numbers: 75.50Lk, 75.40Cx, 75.25+z 
I. INTRODUCTION 

The concept of spin icei has appeared recently as a 
result of a direct analogy between the statistical mechan- 
ics of 0-H bonds in Pauling's model of water ice^'^ and 
those of the localized magnetic moments of some highly 
frustrated magnetai^. These materials, typified by 
Ho2Ti20'ii, Dy2Ti207— and Ho2Sn20'i^, provide almost 
ideal experimental realizations of spin ice, a magnetic 
"sixteen vertex" model''. Their temperature- field phase 
diagrams are exceedingly rich and have revealed several 
new magnetic phases and transitionsiSiSiiSiii. The mag- 
netic dynamics of the spin ices also show some unique 
properties including a crossover from Debye-like spin re- 
laxationiSiiiiJ to quantum tunnelling dynamicsi^^i^ at 
the unusually high temperature of 15 K. The current 
paper builds on earlier worlsiiSiiLi^ to characterize the 
neutron scattering of spin ice materials. In it, we report 
the first neutron scattering images of the spin ice state 

of Dy2Ti207^^ which are compared with theoretical pre- 
diction£2"-2i-22^ 

In Dy2Ti207, as well as in the other spin ice mate- 
rials Ho2Ti207 and IIo2Sn207, the rare earth moments 
('spins') are arranged on a pyrochlore lattice of linked 
tetrahedra (see Fig. P). Strong crystal field anisotropjiS^ 
confines each spin as an Ising-like doublet to the local 



trigonal axis that connects the center of the elementary 
tetrahedron to its vertex. With this constraint and ferro- 
magnetic coupling, the magnetic ground state of a single 
tetrahedron is defined by a "two spins in, two spins out" 
rule, analogous to the Bernal-Fowler ice rule that con- 
trols the proton arrangements in water ice^*. This rule 
minimizes only four of the six near neighbor magnetic 
interactions and so the system is highly frustrated, shar- 
ing with ice the property of zero point entropjiSi^. In 
the spin ice materials, the ferromagnetic near neighbor 
coupling arises from the strong dipole-dipole interaction 
between neighboring spins^. If the dipolar interaction is 
truncated at nearest neighbors one has the "near neigh- 
bor" spin ice modeli*i (note that dipole interactions are 
equivalent to exchange terms for Ising spins). This ap- 
proximation captures the bulk properties of the spin ice 
material9i*^i2Si2S42i and gives a qualitative description of 
the microscopic spin correlations, as determined by neu- 
tron scatteringiiii. 

The subject of this paper, Dy2Ti207, has previously 
been investigated by bulk method o^d^i^^i^^ . In zero field, 
the specific heat resembles a Schottky anomaly, peak- 
ing at 1.24 K and falling to zero at lower temperatures. 
Unlike its IIo2Ti207 analogue, which shows an extra 
peak due to freezing of the nuclear spinsiLSS, Dy2Ti207 
shows no other specific heat contributions down to 0.2 
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FIG. 1: Top: a unit cell of the pyrochlore lattice which con- 
sists of four interpenetrating face centered cubic (fee) Bravais 
sublattices. Two spins related by an fee translation are shown 
with antiferromagnetic correlation. Correlations of this type 
give rise to the zone boundary scattering observed experimen- 
tally. Bottom: the spin ice state on a single tetrahedron (bot- 
tom). A macroscopic spin ice state is constructed by stacking 
such local arrangements on the lattice. 
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TABLE I: Isotopic abundances and absorption cross sections 
{a a) of natural dysprosium and the enriched sample used in 
these experiments^. 



K, where the specific heat becomes exponentially small. 
It has thus proved highly favorable for bulk measure- 
ments, having been used to estimate the zero point en- 
tropy in zero fielc^, and to demonstrate the existence 
of "kagome ice", a field induced state that retains zero 
point entropji2i2£. However, owing to its large neutron 
absorption, Dy2Ti207 has not yet been extensively stud- 
ied by neutron scatteringi2i2£, the technique that was 
crucial in establishing the existence of the spin ice state 
in Ho2Ti207iii and Ho2Sn20A 

To address the problem of large neutron absorption in 
natural dysprosium based materials we have grown an 
isotopically enriched crystal of ^^^Dy2Ti207 specifically 
for use in neutron scattering experiments. The details 
of the sample and the experiments performed on it, are 
given in section FlI Al The experiments were designed to 
examine the nature of the presumed spin ice state that 
is established when the sample is cooled in zero applied 
mag netic field^. The data were modelled using the 
"dipolar" spin ice model that considers further neighbor 
dipole-dipole interactions as well as a small antiferromag- 
netic exchange coupling that can be estimated by fitting 
the susceptibilitji^ or specific heatSS. Details of the dipo- 
lar spin ice model^S, as applied to Dy2Ti207, are given in 
section lilBl of the current paper, while the comparison of 
theory and experiment is given in section UTTl The paper 
is concluded with a discussion of the results, section Hvl 



II. EXPERIMENTAL 
A. Neutron Scattering 

Neutron scattering measurements were carried out 
on PRISMA at the ISIS pulsed neutron source of the 
Rutherford- Appleton Laboratory. A pulsed white neu- 
tron beam is incident on the sample. PRISMA views a 
cold (95 K) methane moderator and a supermirror guide 
system, providing a high flux of long wavelength neu- 
trons allowing measurement at relatively low |Q|. The 
neutrons scattered from the sample are recorded in six- 
teen detectors which cover a range of 16° in scattering 
angle. From the neutron time of flight, the scattering 
angle and the orientation of the sample with respect to 
the incoming beam, the scattered intensity can be deter- 
mined as a function of the position in reciprocal space. 
By rotating the crystal about the vertical axis, inten- 
sity maps covering a large region of the scattering plane 
can be constructed. This method is particularly useful 
for studies of diffuse magnetic scattering. The average 
scattering angle (j) of the sixteen detectors used for the 
measurements is chosen as a compromise between flux 
(at a given wavevector) and background (due to air scat- 
tering). The advantages of measuring the diffuse mag- 
netic neutron scattering on PRISMA have been clearly 
demonstrated using the example of Ho2Ti207, where the 
dipolar spin ice nature of the zero fleld spin correlations 
has been unambiguously establishedii. In that and the 
current experiment the neutron scattering was measured 
in the static approximation'^^. 

Natural dysprosium contains seven isotopes, several of 
which absorb neutrons quite strongly. This can make the 
detection of weaker effects like diffuse scattering difficult 
or impossible. To reduce this, an isotopically enriched 
sample was used. The natural abundances, scattering 
lengths and cross sections for natural dysprosium and 
our sample are given in Table Using the ^^^Dy iso- 
tope the average absorption cross-section CTq is reduced 
by a factor of 4.8 compared to natural Dy. Since absorp- 
tion attenuates scattered intensity by a factor of the form 
exp(— A^CTaA) (where N is the number density of scatter- 
ers per unit volume, aa is the average absorption cross 
section, usually tabulated for 2200 m s~^ neutrons, and A 
is the neutron wavelength scaled to the value used in (Tq) 
this makes a significant difference at long wavelengths. 

A large high quality single crystal was grown using 
an infra-red double mirror image furnace^i. The crys- 
tal has a diameter of « 0.4 cm, length « 1.5 cm, and 
is translucent amber-red in color. It has a cylindrical 
form typical of image furnace grown crystals. The crys- 
tal was varnished into a large copper support to ensure 
good thermal contact. 

Using a '^He sorption refrigerator insert in an Orange 
cryostat, the h, h, I plane of reciprocal space was mapped 
in zero field at 20 K and 0.3 K using the average scatter- 
ing angle (j> = 48°. The same plane was mapped at 0.3 
K and 1.3 K using the average scattering angle 4> = 32°. 
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Using a dilution refrigerator insert in an Oxford Instru- 
ments cryomagnet the h,k,0 and h,h,l planes of recip- 
rocal space were mapped in zero field at 0.05 K. These 
maps are shown in Fig. |21 All parts of the sample en- 
vironment in the beam were made of aluminium (other 
than the copper support), which gives rise to powder lines 
at IQI = 2.70, 3.12 and 4.40 A"!. 

Standard data reduction procedures were applied to 
transform time of flight and angle information to recipro- 
cal space, to correct for absorbtion in the sample, and to 
normalize the data to an absorption corrected vanadium 
run. The latter corrects for the wavelength-dependent 
flux profile and detector efficiency. Despite the five-fold 
reduction in absorption cross section compared to natu- 
ral dysprosium, the data analysis procedure revealed that 
absorption by the Dy2Ti2 07 crystal was still significant 
and needed to be taken into account. The absorption 
correction was made by calculating an attenuation co- 
efficient (A(s,(/))) by which the data could be divided. 
The coefficient for a cylinder of radius R was originally 
derived by Sears^^ and is 

A{s, (j)) = [l + Ahs^ - O.Ss^ cos^{(j)l2)] x exp(-2as), (1) 

where a — (8/3)7r, 6 = (1 — a^)/2, s = fiR and is 
the scattering angle. In this formula, fi, the absorption 
coefficient, depends on both the scattering cross section 
and the wavelength dependent absorption cross section. 



B. Monte Carlo Simulations 

To model the experimental data for Dy2Ti207, we 
use a Hamiltonian appropriate for the (111) Ising py- 
rochlores: 



H =-J 



(2) 



+DRl 



Ca CO 



where Sf — z°-af represents an Ising moment of mag- 
nitude \Sf\ = 1 at fee lattice site i and tetrahedral ba- 
sis position a, with its quantization axis (£°) oriented 
along the local trigonal axis (a member of the (111) set), 
and af — ±1. The nearest neighbor exchange energy 
is Jnn — J/3 and the nearest dipole-dipole strength is 
Dnn = 5D/3, where D = {fio/'lTr)fi^ / R^^ and ^ is the 
moment of the Dy'^^ ion. Rfj is the vector separation 

between moments and S'^. 

The model (|2I), with parameters appropriate to 
Dy2Ti2 07, orders magnetically at 0.18 K^^. The real ma- 
terial does not show this ordering transition, but we test 
the hypothesis that the correlated paramagnetic state 
just above the transition temperature will provide an 
adequate description of the spin ice state. Equation [21 
was simulated on a pyrochlore lattice of 4 x 4 x 4 cubic 



cells {N = 1024 spins) with periodic boundary condi- 
tions. The long range nature of the dipole-dipole inter- 
actions was properly handled by using standard Ewald 
techniques for real space dipoles^S. The model employed 
energy parameters appropriate for Dy2Ti207, J — —3.72 
K and D = 1.41 K^o. 

At temperatures below 1 K or so the system is only 
able to access states within the spin ice manifold. Single 
spin flip dynamics require the introduction of ice rule de- 
fects. This means that the acceptance rate of single spin 
flips becomes extremely small and ability to equilibrate 
the system below 0.3 K becomes dubioua^SiSi. In order 
to improve on approach to equilibrium, dynamics were 
restored by using a loop algorithni2i. To simulate scat- 
tering data the system was placed in the ordered ground- 
state found by Melko et aP^ , at 0.4 K. 5 x 10'^ Monte 
Carlo steps (MCS) per spin and 5 x 10^ loop attempts 
were made to equilibrate the system and bring it to a 
paramagnetic state. Loop attempts were therefore sepa- 
rated by = 1024 MCS. The system was then cooled to 
0.3 K where the same number of single spin flip and loop 
attempt steps were used for equilibration and then data 
collection (see Refs. IT7ll20ll3^ for more details). 

There are 12 symmetry related ordered stuctures and 
to maintain ergodicity each was used as a start point 
for a simulation. I{Q) was calcuated at T = 0.3K from 
200 independent spin configurations derived from each 
starting structure using the formula— 



I{Q) oc 



N 



(iM) 



E 



(3) 



where (...) denotes a thermal average, S""^ is the spin 
component perpendicular to the scattering wave vector 
Q and f{Q) is the magnetic form factor for Dy^+. The 
Monte Carlo results for I{Q) are compared to the exper- 
imental data in the h, h, I and /i, k, planes below. The 
simulations are nominally at 0.3 K, whereas the experi- 
mental data was collected at either 0.3 K (/i, h, I plane) or 
0.05 K (/i, fc, plane). We consider it valid to make these 
comparisons because at temperatures below 0.4 K the 
susceptibilitySi, neutron scattering intensity, I{Q) and 
specific heat^ do not change significantly with temper- 
ature: the system essentially appears to be frozen with 
negligible thermal spin fluctuations. 



III. RESULTS 

In this section, the observed magnetic diffuse scatter- 
ing is presented and compared to Monte Carlo simula- 
tions where possible. The temperature dependence of 
the scattering in the h, h, I plane is discussed first (sec- 
tion llll in order to establish the temperature regimes 
of interest. In the following sections Ijlll Bl and IIII C\i 
the scattering intensity measured at base temperature 
is compared to Monte Carlo simulations. Finally, sec- 
tion is concerned with features in the experimental 
scattering which are not reproduced by the model. 
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A. Temperature Dependence of the Diffuse 
Scattering 

The temperature dependence of the diffuse scattering 
due to the development of ice rule fulfilling spin corre- 
lations was sampled at four temperatures for the h, h, I 
plane (20, 1.3, 0.3 and 0.05 K, see Fig.jSJl. An important 
point that can be obtained from these measurements is 
that the diffuse scattering does evolve with temperature. 
As the sample is cooled, no diffuse scattering character- 
istic of spin ice correlations has built up by 20 K. This 
is entirely expected as the bulk susceptibility shows that 
Dy2Ti2 07 is paramagnetic at this temperatur^. At 1.3 
K strong diffuse scattering has appeared. The form is 
closely similar to that of the near-neighbor spin ice as 
simulated in Ref^3- Between 1.3 and 0.3 K the spin 
correlations evolve and additional features appear in the 
scattering pattern, as discussed below. The scattering at 
0.05 K throughout the h, h, I plane (not illustrated here) 
is not significantly different to that at 0.3 K. 



D. Zone Boundary Features 

One possible explanation for the zone boundary fea- 
tures is that they are due to phonons. However, the fol- 
lowing facts suggest that the intensity is magnetic: (1) 
the scattering appears as temperature is decreased; the 
converse would be expected for a phonon; (2) the exper- 
imentally observed features fall off with increasing |Q|, 
which is not expected for phonons. 

The second point can be demonstrated by comparison 
of the intensity of the zone boundary features with the 
magnetic form factor of Dy"^"*" . If the intensity at equiva- 
lent points on successive zone boundaries is compared it 
would be expected to follow the form factor if the scatter- 
ing process is magnetic. This is found to be the case. On 
the basis of this comparison it seems entirely reasonable 
to attribute the zone boundary features to a magnetic 
correlation. 



IV. DISCUSSION 



B. Scattering in the h,k,Q plane 

For the h, k, plane, data has only been collected at 
0.05 K. The simulation and the experiment seem to agree 
very well indeed. The two scattering maps are shown in 
Fig. 13 A comparison by a cut through the data is shown 
in Fig.0] The Monte Carlo pattern has been fitted to the 
experimental diffuse scattering using a flat background 
and an overall intensity scale parameter. 



C. Scattering in the h,h,l plane 

For the h, h, I plane direct comparison of the experi- 
mental and simulated scattering maps reveals both close 
similarities and significant differences. The two scatter- 
ing patterns are shown in Fig. [21 It is apparent that 
whilst the experimental pattern does contain the features 
of the dipolar spin ice scattering (for example the large 
diffuse feature at 0, 0, 3) it also contains features that 
are are not reproduced by the simulation. These extra 
features are concentrated around the boundaries of the 
Brillouin zones (the position of the zone boundaries are 
illustrated in Fig. |2J|. 

The similarities and differences of the observed and cal- 
culated patterns are more clearly emphasized in Fig. 0] 
Here slices through the data have been made. The agree- 
ment along h,h,0 and 0,0,^ is good. The slice along 
h, h, 2.3 is chosen to highlight the differences and it can 
be seen that where the slice crosses the zone boundary 
the form of the calculated and experimental intensities is 
significantly different. 



The dipolar spin ice model apparently provides a very 
successful model for the zero field spin correlations in 
Dy2Ti2 07. However, some additional spin correlations 
do occur which are not captured by the model. The 
model is very successful in describing the bulk and mi- 
croscopic data for Ho2Ti207ii, and the bulk data for 
Dy2Ti2 07^. One might therefore ask, how this further 
restriction of the dipolar spin ice manifold is not mani- 
fested in the fit to the heat capacity? 

In answer to this question, it should be noted that the 
fit to the heat capacitji^S is not perfect: the experiment 
and simulation differ in detail. The fact that the zone 
boundary features are very diffuse in reciprocal space and 
run all around the zone boundaries, suggests that the as- 
sociated spin correlations are short ranged and isotropic. 
A small difference between the experimental and theo- 
retical heat capacities is therefore understandable since 
the extra spin correlations need not make a major con- 
tribution to the entropy. We therefore conclude that the 
neutron scattering and specific heat datai^iSS are consis- 
tent. 

A clue to the origin of the extra spin correlations is the 
fact that they are much less evident in Ho2Ti207«i2ii&, 
even if they cannot be excluded completely within ex- 
perimental error. According to the model, the principal 
difference between Ho2Ti207 and Dy2Ti207 is the more 
significant antiferromagnetic exchange term in the spin 
Hamiltonian of Dy2Ti207 (J = -3.72 K for Dy2Ti207 
compared to J = —1.56 K for Ho2Ti207). This exchange 
term was invoked to fit bulk measurement3^2*2L^, which 
are not sensitive to the details of a Hamiltonian in the 
way that diffuse neutron scattering is. Hence the extra 
scattering might reveal that extra terms are required in 
the model Hamiltonian: for example, it might be neces- 
sary to include further neighbor exchange interactions. 
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FIG. 2: (Color online) Dy2Ti207: Difltuse scattering in the h,h,l plane measured at 20 K (top left), 1.3 K (top right) and 0.3 
K (bottom left, with zone boundaries illustrated). The sharp, intense spots are nuclear Bragg reflections with no magnetic 
intensity. The map at 20 K was recorded using an average scattering angle oi (j) = 48° which gives stronger arc-shaped artefacts 
for small |Q| values compared to <^ = 32°. The result of the Monte Carlo simulation of the dipolar spin ice model at 0.3 K is 
shown at bottom right. 



FIG. 3: (Color online) Dy2Ti2 07 at 0.05 K: experimental scattering intensity in the h,k,0 plane (left), compared with that 
calculated from Monte Carlo simulations (right). 



Further theoretical work is required to distinguish vari- 
ous possibiUties. 

Zone boundary scattering has recently also been ob- 
served in the h, fc, plane for ZnCr204'^^ and the h, h, I 
plane for ZnFe204— (both frustrated antiferromagnets 
with the spinel structure). We suggest that, while the 
three magnets are quite different, the zone boundary fea- 
tures may also be a generic feature of frustrated mag- 
nets due to the formation of local spin clusters, for ex- 
ample those developing in ZnCr20 4^^'^'^ which might be 
at the origin of the "quantum protectorate" state in- 
voked in that system. As noted by Kamazawa et alm^, 
zone boundary scattering implies correlations between 
the spins related by fee translations (as illustrated in 
Fig. Ql. These spins are third neighbors and interactions 
at this distance will clearly be important in stabilizing a 
hexagonal cluster. Simulations suggest that other clus- 
ters can also be stabilized depending on the details of the 
Hamiltonian'*" . 

In conclusion, we have presented a systematic inves- 
tigation of the diffuse scattering of Dy2Ti2 07 in zero 
magnetic field. We have observed three regimes of be- 
havior. At 20 K the spins are paramagnetic, without 
spin ice correlations, which should help interpret recent 
bulk susceptibility measurements of the dynamical prop- 
gj.^jg^i2ji3j27^ which have been discussed in terms of both 



correlated and uncorrelated dynamical processesi^ti^^i^. 
At 1.3 K the spin correlations are similar to those of 
the near neighbor spin ice model, suggesting that further 
neighbor dipole interactions are not strongly relevant in 
this temperature range. At 0.3 K the scattering is modi- 
fied to correspond quite closely to the frozen state of the 
dipolar spin ice model, with effective infinite range in- 
teractions. However, in this temperature range we have 
observed additional scattering features around the zone 
boundary which for the time being remain unexplained. 
Further calculations are in progress to clarify their origin. 
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